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Develop  an  active  sonar  (Lorad)  for  the  detection  of  sub- 
marines in  deep  water  at  ranges  of  100  miles  or  more, 
and  in  shallow  water  in  excess  of  10  miles.  Specifically, 
design,  install,  and  evaluate  an  improved  receiving  sys- 
tem for  Lorad  Mark  III. 


RESULTS 


A receiving  system  designated  DIMUS  (Digital  MUltibeam 
Steering)  was  developed.  Tests  showed  that  this  system, 
using  five  samples  per  cycle,  has  satisfactory  beam-form- 
ing capability  and  its  bearing  precision  is  equal  to  that 
obtainable  with  other  methods,  such  as  delay  lines. 


RECOMMENDATION 


Use  DIMUS  in  place  of  delay  lines  in  future  Lorad  design, 
provided  that  the  simpler,  passive  phase  compensator 
cannot  be  employed. 
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SYMBOLS 


n = Element  number  from  1 to  60 
d = Spacing  between  elements 

9 = Angle  of  arrival  measured  from  array  normal 

9C  = Angle  of  compensation 
AT  = Time  delay 
c = Sound  velocity 

Bw  = Number  of  bits  assigned  to  water  path 
Bs  = Number  of  bits  in  shift  register 
k = Integer 

A = Positive- going  point  on  square  wave 

C = Array  center 

R = Amplitude  of  beam  output 

P - Polarity  of  1500- c/s  square  wave 

Bn  - Bit  difference  relative  to  array  center 

F = Clock  frequency 

N - Number  of  shift  register  bits 
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INTRODUCTION 


Beam  formation  in  the  past  has  usually  been  accomplished 
by  means  of  delay  lines.  In  a steered  sonar  system,  differ- 
ent sets  of  adjacent  hydrophones  are,  in  turn,  connected  to 
the  same  delay  line,  thus  producing  maximum  response 
throughout  360°  in  azimuth,  but  in  only  one  direction  at  a 
time.  Hence,  the  manual  search  rate  is  rather  low.  If  a 
mechanical  scanning  switch  is  employed,  time- sharing  and 
loss  of  information  result. 

In  a multibeam  sonar,  a delay  line  can  be  used  for  each 
desired  beam,  and  each  hydrophone  output  is  connected  to 
the  proper  tap  on  each  delay  line.  Since  all  outputs  are 
available  to  the  operator  simultaneously  there  is  no  search- 
rate  problem,  and  all  available  information  is  fully  used. 

An  example  is  an  installation  made  for  Lorad  aboard 
USS  BAY  A.1 

Mechanical  scanning  switches  introduce  problems  of 
switching  noise  and  maintenance,  while  batteries  of  delay 
lines  take  up  scarce  shipboard  floor  space. 

Another  approach  to  multibeam  formation  is  phase  compensa- 
tion, which  involves  the  use  of  transformers  with  pairs  of 
secondary  windings  equaling  the  number  of  beams  to  be 
formed.  Such  a phase  compensator  was  developed  and 
installed  for  Lorad  station- keeping.  2 Phase  compensators 
are  primarily  intended  for  single- frequency  or  narrow- 
band  operation  but  their  applicability  to  broad-band  opera- 
tion has  also  been  demonstrated.  ^ 


^Navy  Electronics  Laboratory  Report  968,  Lorad  Mark  III 

Receiving  System- Bearing  Calibration  and  Self- Noise,  by 

Cl  T.— Krieger,  CONFIDENTIAL,  'll  May  1960 
o 

Navy  Electronics  Laboratory  Report  757,  Simultaneous 
Multibeam  Phase  Compensation,  Part  VI:  Phase  Compensa- 
tor for  the  Lorad  Station- Keeping  Receiver  System,  by  C.  J. 
Krieger  and  R P.  Kempff,  CONFIDENTIAL,  7 May  1957 

3 

Navy  Electronics  Laboratory  Report  1009,  Simultaneous 
Multibeam  Phase  Compensation:  Part  X Broadband 
Application  by  F.  KT  Abbott  8 December  1960 
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A third  approach  to  multibeam  formation  is  DIMUS  (Digital 
MUltibeam  Steering).  The  clipped  signals  are  fed  into 
shift  registers,  and  advanced  by  means  of  clock  pulses. 

By  connecting  the  proper  taps  on  the  shift  registers,  any 
number  of  desired  beams  can  be  obtained  simultaneously 
(fig.  1). 

The  study  to  be  reported  here  includes  the  design  of  DIMUS, 
laboratory  evaluation  and  sea  trials,  and  a comparison  of 
its  performance  with  that  of  delay  lines. 


^ Anderson,  V.  C.  , "Digital  Array  Phasing,  " Acoustical 
Society  of  America.  Journal,  v.  32,  p.867-870,  July  1960 


BEAM  BEARING  BUS 


Figure  1.  Schematic  diagram  of  Lorad  DIMUS. 
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DETAILS  OF  DIM  US  DESIGN 


EXPRESSIONS  FOR  THE  DELAYS  IN  THE  WATER  PATH 
AND  THE  SHIFT  REGISTER 


Delay  in  the  Water  Path  Referred  to  the  Array  Center 


It  may  be  seen  from  figure  2 that  the  path  difference  at  the 
nth  element,  referred  to  the  array  center,  is  (30.  5 - n ) 
d • sin  S,  where 

n is  the  element  number  from  1 to  60 

d the  spacing  between  elements 

Q the  angle  of  arrival  measured  from  the 
array  normal 

The  time  delay  AT  at  the  n n element,  still  referred  to  the 
array  center,  is  then 

(30.  5 - n)  d • sin  Q 
c 

and  if  d = 20  inches  or  1-2/3  feet,  and  sound  velocity 
c = 4900  ft/sec,  AT  = (30.  5 - n)  sin  0 • 340  • 10'°  sec. 


Figure  2.  Path  difference  of  plane  wave  referred  to  array 
center  at  nth  element  of  a linear  array. 
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Let  us  now  assume  a signal  frequency  of  1500  c/s  and 
assume  that  the  information  is  advanced  along  the  shift 
register  at  the  rate  of  7 500  clock  pulses  per  second;  then 
the  corresponding  number  of  bits  in  the  water  path  is: 

(30.  5 - n)  sin  9 • 340  • 7500  • 10‘6  = 2.  55  sin  9 (30.  5 - n) 

where  one  bit  is  equivalent  to  one  flip-flop  in  the  shift 
register.  To  eliminate  negative  delays,  which  cannot  be 
realized  with  this  method,  we  add  the  maximum  number 
of  bits  expected  at  either  array  end  for  0 = ±60°,  namely 

29.  5 • sin  60°  • 340  • 7500  • 10-6  = 65.  2 

Hence,  the  number  of  bits  B assigned  to  the  water  path  is 

Bw  = 65.  2 + 2.55  • sin  9 • (30.  5 - n) 


Delay  in  the  Shift  Register  Referred  to  the  Array  Center 


In  order  to  refer  the  arrival  time  from  all  directions  to 
the  array  center,  and  since  the  maximum  angles  of  com- 
pensation are  ±60°,  the  number  of  bits  in  the  shift  register 
Bs  are  chosen  so  that  Bw  + Bs  = 130.  4 when  0 = 0C  for  any 
hydrophone  n.  (ec  is  the  angle  of  compensation.  ) This 
demands  that 

B - 65.2  + 2. 55  • sin  9. • (n  - 30. 5) 

S b 

One  stage  (bit)  was  added  at  the  input  of  each  shift  register 
in  order  to  provide  a tap  for  zero  delay. 


Combined  Delay  in  the  Water  Path  and  in  the  Shift  Register 


The  total  number  of  bits  is  then 

Bw  + Ba  = 1 30.  4 + 2.  55  (sin  9C  - sin  e)  • (n  - 30.  5) 

In  choosing  shift  register  taps,  fractions  are  rounded  up 
or  down,  and  the  actual  number  of  taps  may  differ  from 
the  calculated  number  by  ±1. 
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DIRECTIVITY  PATTERN  WITH  SINGLE  FREQUENCY 
SIGNAL 


Maximum  output  is  obtained  if  all  shift  register  taps  used 
in  beam  formation  are  in  the  same  state,  say  "l,  ' at  one 
time.  A given  tap  delivers  information  either  "l"  or  "0,  " 
7500  times  per  second.  If  the  shift  register  inputs  are 
1500- c/s  square  waves,  then  either  three  "l's"  or  two 
"0's,  " or  the  reverse,  will  be  delivered  at  any  tap  during 
one  cycle,  i.e.,  in  i/1500  second.  The  amplitude  variation 
in  time  of  one  beam,  holding  ec  and  d constant,  can  be 
obtained  analytically  by  adding  the  outputs  of  all  taps 
involved  in  forming  this  beam,  at  selected  instants  during 
one  cycle. 

There  is  a simple  way  of  finding  the  elements  n where  the 
difference  of  bits  Bw  + B s relative  to  the  array  center  has 
reached  ±1/2,  ±3/2,  ±5/2,  etc.,  and  where,  therefore,  a 
whole  bit  must  be  added  or  subtracted. 

Set 

2.55  (sin  ec  - sin  q)  (n  - 30.  5) 

= ±1/2,  ±3/2 ~~2~’  *s  1-2.  . . 

Solving, 

_ c . 2 k-  \ 

71  — 30.  5 — p-"  ' -C  r\  f . " ^ 

5.  10(sin  Qc  - sin  0) 

For  6 = 6C,  no  values  of  n between  1 and  60  exist,  and 
therefore  all  taps  deliver  the  same  information. 


Now  let 

0C  = 60°  and  e = 54°,  with  sin  9C  - sin  Q - .05701, 

» = 30-  5 ± 2910  = 30.  5 ± 3.  45  or  27.  05  and  33.  95 

and 

n = 30.  5 ± 29iq  = 20.  2 and  40.  7,  and  so  forth. 
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If  we  refer  all  bits  to  the  array  center,  then  the  distribution 
of  bits  shown  in  figure  3 results. 


BIT  DIFFERENCES  REFERRED  TO  ARRAY  CENTER 
-4  I -.3  I -2  I -I  I 0 1+1  1+2  1+3  I +4 


1 6,7  13,14  20,21  27,28  33,  34  40,41  47,48  54,55  60 

677  767776 

NUMBER  OF  ELEMENTS  WITH  A SPECIFIED  BIT  DIFFERENCE 

Figure  3.  Distribution  of  bit  differences  referred  to  array 
center,  for  oc  - 60°  and  y = 54°. 


The  algebraic  sign  in  the  left  half  of  the  array  is  taken  as 
negative  since  at  q = 54°  the  water  path  there  is  shorter 
then  at  uc  = 60°.  These  bit  differences,  it  should  be  realized, 
are  caused  by  the  fact  that  the  plane  wave  arrives  from  a 
direction  other  than  the  direction  of  compensation,  and  that 
the  plane  wave  is  sampled  five  times  per  cycle.  This  results 
in  either  wo  "i"  states  and  three  0"  states  in  one  cycle  of 
the  signal  frequency,  or  the  reverse. 

The  total  output,  at  any  time,  will  depend  on  the  relative 
position  of  the  advancing  square  wave  which  carries  the 
information.  Let  us  examine  the  situation  when  the  positive- 
going point  A on  the  square  wave  has  just  passed  the  array 
center  C,  where  the  bit  difference  is  zero  (fig.  4).  At  this 
moment  there  will  be  33  "l's";  two- tenths  of  a cycle  later, 
in  the  dotted  position,  there  will  be  32  "l's.  " Continuing 
this  process  we  arrive  at  the  following  amplitude  variation: 

Phase  of  square  wave  at  c 
counted  from  positive- going 

point  A in  fractions  of  Amplitude  of 

signal  period beam  output  ( R ) 


0 

33 

0.2 

32 

0.  4 

33 

0.  6 

41 

0.  8 

41 

1.0 

33 
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angles  of  compensation  and  several  angles  of  arrival,  the 
results  shown  in  figures  5 and  6 were  obtained.  The  result- 
ing curves  contain  the  1500- c/s  signal  frequency  in  which 
we  have  the  greatest  interest,  in  addition  to  a constant  and 
some  higher  frequencies.  The  former  can  be  isolated  by 
means  of  a band- pass  filter  or  found,  approximately  at  least, 
by  a graphical  analysis  as  shown  in  the  figures. 


The  maximum  amplitudes  of  the  signal  are  a measure  of  the 
response  of  the  system  to  a signal  arriving  from  the  chosen 
direction,  i.e.,  the  directivity  function.  Two  of  these  have 
been  plotted  in  figures  7 and  8;  the  values  are  shown  in 
solid  lines.  There  is  good  agreement  with  the  values 
obtained  by  the  usual  analog  method,  employing  delay  lines. 
The  corresponding  values  are  plotted  in  dotted  lines. 


n 

If  ten  samples  instead  of  five  were  taken  in  one  cycle,  the 
clock  frequency  and  the  number  of  bits  per  element  would 
be  doubled  for  the  same  bearing  angle  coverage.  The 
resulting  responses  are  shown  in  figures  9 and  10;  the 
corresponding  points  are  shown  as  dashed  lines  in  figures 
7 and  8.  The  results  do  not  differ  appreciably  from  those 
obtained  by  taking  five  samples;  the  latter  choice  is  there- 
fore justified. 


Analytical  Method 


4 

A method  has  been  devised  which  can  be  adapted  for  use 
with  a computer,  if  desired.  In  brief,  first  the  bit  difference 
relative  to  the  array  center  is  calculated  for  the  taps  which 
are  used  in  beam  formation.  Then  the  upward  crossing  of 
the  square  wave  representing  the  information  is  lined  up 
with  the  array  center,  and  the  polarity  of  the  square  wave 
is  noted  at  all  taps  involved  in  beam  formation.  The  sum 
of  all  positive  polarities  is  a measure  of  the  response. 

Next,  the  square  wave  is  advanced  one  bit,  or  2 it/ 5 radians 
in  phase  at  a time,  and  the  process  is  repeated  until  five 
samples  have  been  obtained  for  each  direction  of  compensa- 
tion and  angle  of  arrival. 


Bit  Differences  Relative  to  Array  Center 

From  an  earlier  equation  we  have 

Bw  + Bg  = 1 20.  4 + 2 . 55  (sin  0C  - sin  o)  {n  - 30.  5) 
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Figure  7.  DIMUS  directivity  pattern  as  a function 
of  6;  0C  = 0°;  five  and  ten  samples  per  period, 
also  analog  process;  1500-c/s  signal  frequency. 
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Figure  9.  Amplitude  variation  for  0C  = 0°  a 
various  values  of  0;  ten  samples  per  period; 
1500-c/s  signal  frequency. 
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For  our  purpose  we  are  interested  only  in  the  bit  differences 
relative  to  the  array  center,  where  Bw  + Ba  = 130.  4.  Let 
us  then  define  the  bit  differences  relative  to  the  array  center 
as 

B„  = 2.55  (sin  er  - sin  g)  • (n  - 30.  5) 

71 

where,  as  before,  ec  and  g are  the  angles  of  compensation 
and  arrival,  respectively,  and  n is  the  element  number  from 
1 to  60.  Since  Bn  is  quantized,  it  must  be  rounded  up  or 
down  to  the  nearest  integer. 


Sampling  of  Polarity  of  Square  Wave 


' 


The  1500-c/s  square  wave  which  carries  the  information 
advances  down  each  shift  register  at  the  rate  of  one  cycle 
per  five  clock  pulses.  Hence  the  polarity  of  the  square 
wave  can  be  sampled  five  times  per  cycle.  This  occurs 
simultaneously  in  each  of  the  60  shift  registers. 

Let  us  examine  a plane  sound  wave  arriving  at  the  same 
time  at  all  the  elements  of  an  uncompensated  array.  The 
polarities  at  the  corresponding  taps  will  be  the  same,  and 
the  signals  will  enhance  each  other. 

If  now  the  sound  wave  arrives  at  hydrophone  No.  1 before 
it  arrives  at  the  array  center,  the  water  path  is  shorter 
there  and  the  information  appearing  in  shift  register  No.  1 
is  earlier  than  that  at  the  same  tap  at  the  array  center. 

Let  us  designate  this  as  a negative  bit  number. 

If,  for  example,  the  square  wave  is  sampled  just  after  the 
upward  crossing  at  the  array  center,  then  a positive  polarity 
is  found.  (See  diagram  below.  ) If, at  the  same  moment,  it 
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is  sampled  at  hydrophone  No.  1 where  the  bit  difference 
may  be  -4,  then  the  polarity  would  be  0.  At  hydrophone 
No.  60,  on  the  other  hand,  where  the  bit  difference  may 
be  +4,  the  polarity  would  be  1.  The  sum  of  the  last  two 
outputs  would  thus  be  1.  The  58  other  outputs  would  have 
to  be  added  to  obtain  the  total  array  response. 

We  can  summarize  our  findings  in  the  form  of  a table: 


+5 

+4 

+3 

+2 

+1 

0 

-1 

-2 

-3 

-4 

-5 

1 

1 

1 

0 

0 

1 

1 

1 

0 

0 

1 

One- fifth  of  a cycle,  or  one  pulse,  later,  more  recent 
information  appears  at  all  taps.  This  can  be  expressed 
by  stating  that  the  bit  differences  Bn  have  increased  by 
-1  (or  decreased  by  1).  It  can  also  be  interpreted  as  a 
motion  of  the  square  wave  to  the  left,  or  down  the  shift 
register: 


The  System  Response 


Proceeding  as  before,  the  responses  in  various  directions 
are  found. 

Table  1 illustrates  the  procedure,  and  results  for  Qc  = 60° 
and  0 = 54°  are  shown.  The  column  headings  are 

n = element  number 

Bn  = bit  difference  relative  to  array  center 
Bc  = bit  number  at  array  center 
P = polarity  of  1500  c/s  square  wave 
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Bn  = 0.  1455  (n  - 30.  5)  and  the  system  response  is  if  = 

n-  1 

The  system  response  if  ranges  from  32  to  41  and  is,  as 
might  be  expected, identical  with  the  results  obtained  with 
the  nomographic  method  (fig.  6). 

In  general,  the  nomographic  method  Will  be  preferable  if 
only  a small  number  of  points  on  the  directivity  pattern 
are  wanted;  on  the  other  hand,  the  analytical  method  when 
adapted  for  computer  use  would  be  chosen  if  a large  number 
of  points  are  wanted. 

In  an  actual  system,  it  may  be  expected  that  fluctuations  in 
the  relative  phases  at  the  array  elements,  finite  rise  time, 
and  other  departures  from  the  underlying  assumptions, 
will  result  in  a statistical  rather  than  strictly  algebraic 
summing  of  the  sixty  tap  outputs. 


DIRECTIVITY  PATTERN  WITH  BROADBAND  SIGNAL 


The  preceding  discussion  has  shown  that  the  directivity 
pattern  obtained  from  a single  frequency,  using  shift 
registers,  is  practically  identical  with  that  obtained  with 
delay  lines.  If  another  single  frequency  differing  from  the 
design  frequency  is  used,  then  the  number  of  samples 
taken  per  cycle  will  also  differ,  but  the  previous  discussion 
and  figures  9 and  10  indicate  that  the  direction  of  maximum 
response  is  not  altered.  It  is  further  concluded  that  this 
statement  holds  also  for  a broadband  signal,  which  is 
composed  of  numerous  single  frequencies.  However, 
changes  in  minor  lobes  and  nulls  must  be  expected. 

To  illustrate  the  results  that  may  be  expected,  a signal 
composed  of  three  frequencies,  1450,  1500,  and  1550  c/s, 
and  a sampling  frequency  of  7500  c/s  were  chosen,  and 
the  resulting  directivity  pattern  was  obtained  with  the 
nomographic  window  method.  Table  2 shows  the  responses 
at  different  directions  of  arrival  for  the  three  frequencies. 
The  amplitudes  R shown  in  the  last  column  were  obtained 
from  graphs  or  by  inspection  of  the  tabulated  values;  the 
response  of  36  at  60°,  however,  was  obtained  by  a Fourier 
analysis.  The  responses  at  the  three  frequencies  were 
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TABLE  2.  SYSTEM  RESPONSE  R FOR  Qc  = 60°,  DIMUS 
CLOCK  FREQUENCY  7500  C/S. 

SIGNAL  FREQUENCIES  1450,  1500,  AND  1550  C/S 


PHASE 


9 

-.8 

-.6 

-•4  1 

-.2 

o 

| +.2  | 

+ . 4 

+ .6  j 

i +.  8 i 

R 

14 

50  C/S 

60° 

60 

60 

00 

0 

60 

60 

60 

0 

0 

+ 36* 

58° 

60 

41 

19 

19 

41 

60 

41 

19 

0 

+ 30 

57° 

44 

37 

! 31 

31 

37 

44 

37 

23 

16 

+ 14 

56°25' 

36 

36 

36 

36 

36 

36 

36 

! 24 

24 

+ 6 

56° 

32 

35 

38 

39 

35 

32 

31 

I 29 

28 

-5.  5 

55° 

26 

34 

43 

43 

34 

25 

26 

34 

35 

-9 

54° 

26 

33 

41 

1 41 

33 

1 26  j 

25 

34 

1 34 

-8 

1500  C/S 


60° 

60 

60 

0 

0 

60 

60 

60 

0 

0 

+ 36* 

58° 

60 

41 

19 

19 

41 

60 

41 

19 

19 

+20.  5 

57° 

44 

37 

31 

31 

37 

44 

37 

31 

31 

+ 6.5 

56°25' 

36 

36 

36 

36 

36 

36 

36 

36 

36 

0 

56° 

32 

35 

39 

39 

35 

32 

35 

39 

39 

-3.  5 

55° 

26 

34 

43 

43 

34 

26 

34 

43 

43 

-8.  5 

54° 

32 

33 

41 

41 

33 

32 

33 

41 

41 

-4.  5 

1550  C/S 


60° 

; 60 

60 

1 0 

0 

60 

60 

60 

1 

0 

0 

+ 36* 

58° 

41 

41 

19 

19 

41 

60 

41 

19 

19 

+20.  5 

57° 

34 

29 

31 

31 

37 

44 

37 

31 

31 

+ 7.  5 

56°25' 

24 

24 

36 

36 

36 

36 

36 

36 

36 

-6 

56° 

21 

25 

35 

39 

35 

32 

35 

39 

39 

-9 

55° 

18 

25 

35 

42 

34 

26 

34 

43 

43 

-12.  5 

54° 

25 

26 

34  j 

35 

33 

32 

33  j 

41 

41 

-8 

*By  Fourier  analysis 
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combined  assuming  a random  phase  relationship;  the  resulting 
directivity  pattern  obtained  is  shown  in  table  3.  The  pattern 
is  presented  in  figure  1 1 as  a broken  curve.  As  predicted, 
maximum  response  is  still  at  60°;  the  minor  lobe  is  about 
- 1 1 db  and  the  depth  of  the  first  null  is  about  - 17  db. 


j 


TABLE  3.  DIRECTIVITY  PATTERN  R FOR  6C  = 60°,  DIMUS 
CLOCK  FREQUENCY  7500  C/S. 

SIGNAL  FREQUENCIES  1450,  1500,  1550  C/S 


Q 

2 

R 

| 2 R2 

R 

Rel. 

db 

1450 

C/S 

1500 

C/S 

1550  1 
C/S 

db 

Vz  r£ 

60° 

1300 

1300 

1300 

3900 

62.5 

1 

35.  9 

0 

58° 

900 

420 

420 

1740 

41. 8 

32. 4 

-3.  5 

57° 

196 

42 

56 

294 

17.2 

24.  7 

-11.2 

56°25' 

36 

0 

36 

72 

8.  5 

18.  6 

- 17.  3 

56° 

30 

12 

81 

123 

11.  1 

20.  9 

- 15.  0 

55° 

81 

72  j 

156 

309 

17.  5 

24.  9 

-11.0 

54° 

64 

20 

1 

64 

148 

12.2 

21.7 

-14.  2 

For  comparison  purposes,  the  directivity  pattern  resulting 
from  a delay  line  was  then  computed  from  the  expression 

sin  60  — (sin  6 „ - sin  6) 


*o  sin  — (sin  Qn  - sin  6) 

\ c 

where  Qc  = 60°  and  d - 5/3  foot.  The  same  three  frequencies 
and  random  phase  relationship  were  used,  as  above.  The 
resulting  directivity  pattern  is  shown  in  figure  11  as  a solid 
curve.  The  general  character  of  the  two  patterns  is  quite 


26 


CONFIDENTIAL 


* Mi 


CONFIDENTIAL 


BEARING  (°> 


Figure  11.  Comparative  directivity  patterns  for  DIMUS 
and  delay  line  as  a function  of  q;  oc  ~ 60°;  signal  frequencies 
1450,  1500,  and  1550  c/s. 


similar,  but  the  delay  line  pattern  has  a lower  minor  lobe 
level  and  a deeper  null.  It  should  be  borne  in  mind  that  in 
the  DIMUS  analysis  a signal  duration  of  less  than  2 cycles 
was  used;  for  longer  signals,  a time  dependency  of  the  pattern 
may  be  expected,  but  an  averaged  pattern  will  result. 


SHIFT  REGISTERS  AND  BEARING  BUS  ASSEMBLY 


No  attempts  are  made  in  this  report  to  present  the  details  of 
equipment  design.  While  the  preceding  discussion  is  based 
on  well  known  conventional  shift  registers,  it  applies  equally 
well  to  the  shift  register  actually  used  in  this  application. 

The  operation  of  these  flip-flops  is  unique  because  informa- 
tion is  advanced  by  means  of  an  overlapping  two- phase  clock 
which  turns  the  collector  supply  on  and  off.  This  results  in 
a flip-flop  output  which  is  the  amplitude- modulated  clock 
pulse  sequence  (fig.  12). 

Although  the  waveforms  at  the  shift  register  taps  are  different 
from  those  in  the  conventional  circuits,  the  methods  of  cal- 
culating directivity  patterns  given  earlier  are  still  applicable. 
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"A"  CLOCk'1 


L OVERLAPPING 
f TWO- PHASE  CLOCK 


"B"  CLOCK J 


INPUT  SIGNAL 
1 500-  C 1 S CLIPPED 
SINE  WAVE 


OUTPUT  OF  A FLIP-  FLOP 
IN  SHIFT  REGISTER 
DRIVEN  BY  "B"  CLOCK 


Figure  12.  Waveform  of  flip-flop  using  overlapping  two- 
phase  clock. 


Through  the  cooperation  of  the  Marine  Physical  Laboratory 
of  the  University  of  California,  the  required  shift- register 
boards  were  procured.  The  taps  to  be  connected  for  the 
formation  of  the  121  beams  were  determined  by  the  methods 
described  previously  (p.  9 ) and  tabulated.  The  wiring  of 
the  supporting  rails  was  accomplished  by  an  outside 
contractor;  final  assembly  was  performed  at  NEL.  Fig- 
ure 13  shows  the  DIMUS  cabinet  installed  aboard  USS 
BAY  A. 

Detailed  wiring  diagrams  are  not  included  in  this  report, 
but  are  available  upon  request.  Address  U.  S.  Navy 
Electronics  Laboratory,  San  Diego  52,  California,  atten- 
tion Code  2 320. 
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Figure  13A.  DIM  US  cabinet  installation  aboard 
USS  RAYA,  with  front  panel  in  place. 


CONFIDENTIAL 


j <t 


CONFIDENTIAL 


LABORATORY  EVALUATION 

RESPONSE  TO  SIMULATED  PLANE  WAVES 

Upon  completion  of  the  121- beam  DIMUS  matrix,  input 
signals  simulating  arriving  plane  waves  were  applied  to 
the  inputs,  after  clipping.  The  matrix  output  was  passed 
through  a 1250-  to  1750-c/s  band-pass  filter.  Three  test 
results  are  described  here: 

A.  An  input  signal  of  1500  c/s  was  applied  in  phase  to 
56  operative  channels  (of  the  60  design  channels).  This 
test  simulates  a target  at  relative  bearing  000°. 

Figure  14  shows  the  response  of  the  121  output  channels 
or  bearing  buses,  numbered  from  1 to  121  to  correspond 
to  relative  bearings  from  300°  to  060°.  It  should  be 
emphasized  that  this  is  not  a directivity  pattern;  rather 
it  depicts  the  simultaneous  outputs  of  the  121  channels, 
in  response  to  a source  at  000°.  Maximum  response 
does  occur  in  channel  61,  i.e.,  at  relative  bearing  000°. 
The  response  of  the  adjacent  channels  is  about  4 db  lower, 
as  expected  from  their  beam  width  of  about  2°;  beam 
separation  is  1°.  All  other  channels  respond  substantially 
less,  in  keeping  with  the  design  objectives. 

B.  This  test  was  the  same  as  test  A,  except  that  the  input 
signal  was  random  noise  in  a 1200-  to  1700-c/s  band,  in 
phase. 

Figure  15  shows  that  the  response  at  the  bearing  of  the 
source  (000°)  and  of  the  adjacent  channels  is  quite  similar 
to  that  produced  in  test  A.  The  response  of  the  remaining 
channels  is  more  uniform  and  lower,  without  marked 
spikes.  This  statistical  smoothing  effect  is  to  be  expected 
with  broad- band  inputs. 

C.  An  input  signal  of  random  noise  in  a 1200-  to  1700-c/s 
band  was  applied  to  channels  1 to  50  with  required  phase 
differences,  by  means  of  a delay  line,  to  simulate  a plane 
wave  arriving  from  relative  bearing  311.5°.  (This  bearing 
resulted  from  the  use  of  an  available  delay  line  designed 
for  different  element  spacing.  ) 

Figure  16  shows  that  maximum  responses  occur  at  relative 
bearings  311°  and  312°;  hence  the  relative  bearing  of  the 
source  is  311.5°,  as  expected.  The  responses  of  the 
remaining  channels  are  uniform  and  low,  as  in  test  B. 
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A clock  frequency  of  7.5  kc/s  ± 1 c/s  was  maintained  by 
monitoring  with  an  electronic  counter,  since  the  indicated 
bearing  is  a function  of  the  clock  frequency. 


EFFECT  OF  CLOCK  FREQUENCY  ON  INDICATED 
BEARING 


The  number  of  shift  register  bits  N needed  between  adjacent 

elements  is  — ~ - F where  p is  the  clock  frequency,  and 

d,  c,  e have  their  usual  meaning.  Consider  a plane  wave 
arriving  from  300°  at  the  extreme  elements,  numbers  1 
and  60  of  the  array.  The  difference  in  bits  required  is 
59  d 

N = sin  Q • F.  If  F is  increased  to  F',  then  the  informa- 

C C r q , 

tion  in  channel  1 has  advanced  through  N'  = — — sin  d • F ' 

c 

units  by  the  time  the  same  plane  wave  arrives  at  element 
60.  Hence  phasing  is  accomplished  at  a bearing  bus  bar 
which  was  designed  for  a bearing  e at  the  normal  clock 
frequency  F.  Thus, 


59  d . 5 9 d . 

sin  0 • f 1 = sin  q.  • f 

c c c 

and  the  indicated  bearing  ec  is  defined  by 


sin  9.  = sin  e • — 
c F 

i.  e. , farther  from  the  array  normal  than  the  correct 
bearing  9. 


The  following  values  were  obtained  experimentally  when 
the  source  was  simulated  at  0 = 311.  5 with  a design  clock 
frequency  F = 7.5  kc/s: 


(kc/s) 

Measured 

Sc 

Calculated 

0C 

7.4 

312.5 

312.4 

7.  45 

312 

311. 9 

7.  5 

311.  5 

311.  5 

7.  6 

310.  5 

310.  6 
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The  agreement  is  very  satisfactory.  Clearly,  it  is  essen- 
tial to  use  a very  stable  and  controllable  clock  frequency. 

At  the  same  time,  the  relation  between  clock  frequency 
and  indicated  bearing  can  be  utilized  for  beam  steering. 

From 


N = 1L1  . sin  6 • F 
c c 


we  find 


ddc  = - tan  ec  • f- 
provided  that  c is  constant. 


Hence  a small  change  in  clock  frequency  will  result  in  a 
predictable  change  in  the  indicated  bearing. 


CLOCK  FREQUENCIES  FOR  DIFFERENT  VELOCITIES 
OF  SOUND 


It  can  also  be  shown  that  the  indicated  bearing  Qc  will 
change  if  the  sound  velocity  c differs  from  the  design 
value  while  the  design  clock  frequency  F is  constant.  In 
order  to  maintain  a correct  dc , F must  be  changed.  From 
the  equation 


N = 5 — sin  q 
c 


F 


it  follows  that 


F 1 = 


c'F 

c 


where  F'  and  c'  are  the  new  values  of  clock  frequency  and 
sound  velocity,  respectively,  and  F and  c are  the  corre- 
sponding design  values. 
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BEARING  CALIBRATIONS 


Bearing  calibrations  were  made  at  sea  on  26  and  27  January 
1961.  USS  REXBURG  (EPCER  855).  while  lying  to  in  700 
fathoms,  lowered  a 2 4- inch- diameter  barium  titanate  trans- 
ducer 50  feet  and  transmitted  0.  5- second  pulses  spaced 
2.0  seconds  at  1500  c/s,  or  random  noise  in  a band  from 
1250  c/s  to  1750  c/s.j  USS  BAYA  (AGSS  318),  on  which  a 
100- foot  linear  array  and  the  DIMUS  system  were  installed 
described  tight  counterclockwise  circles  at  periscope  depth. 
The  average  range  was  approximately  2000  yards;  range 
measurements  were  made  three  times  during  each  run. 
Periscope  bearings  of  the  source  were  obtained  at  1°  inter- 
vals and  marked  on  two  paper  recorders,  together  with 
the  outputs  of  ten  bearing  channels  (also  called  beams). 
Consequently,  both  the  acoustical  and  periscope  bearings 
were  known  simultaneously. 

Five  runs  were  made  with  the  DIMUS  equipment,  and  500 
bearing  readings  were  taken,  i.e.  , an  average  of  about 
four  for  each  of  the  121  beams. 

Also,  three  runs  were  made  with  the  delay  lines.  ^ In  one 
run,  the  outputs  of  the  DIMUS  and  the  delay  lines  were 
recorded  simultaneously,  for  a direct  comparison. 

The  results  of  the  five  DIMUS  runs  are  shown  in  figure  17; 
each  run  is  represented  by  symbols  and  colors. 

The  periscope  readings  have  not  been  smoothed,  although 
their  irregular  spacing  in  time  suggests  a fair  amount  of 
scatter.  A smoothing  procedure  would  involve  an  unjustified 
expenditure  of  time  and  effort. 

The  acoustical  bearings  have  been  corrected  for  the  two 
known  effects,  namely  parallax  and  velocity  of  sound.  The 
array  is  located  about  100  feet  forward  of  the  periscope 
and  the  parallax  is  of  the  order  of  about  0.  5°,  depending 
on  both  range  and  bearing.  A sound  velocity  of  4932  ft/sec 
was  obtained  from  bathythermograms  taken  several  times 
during  the  test.  Corrections  to  the  design  velocity  of 
4900  ft /sec  are  of  about  the  same  magnitude  as  the  parallax 
corrections,  but  have  the  opposite  algebraic  sign. 

Figure  17  shows  that  the  errors  are  independent  of  bearing, 
and  that  the  majority  carry  a negative  algebraic  sign.  The 
extremes  are  +1.4°  and  -1.4°. 
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The  distribution  of  the  errors  is  shown  in  figure  18.  As 
stated  before,  corrections  for  parallax  and  velocity  of 
sound  had  been  made.  The  mean  x = 0.  18°  indicates  that 
periscope  readings  could  be  increased  by  +0.  18°  or  11' 
to  make  the  periscope  zero  point  agree  with  the  average 
acoustical  zero  point.  Conversely,  all  acoustical  bearings 
could  be  decreased  by  the  same  amount.  The  standard 
deviation  is  +0.45°;  this  includes  errors  in  reading  peri- 
scope bearing,  and  in  locating  maximum  response  on  the 
paper  records;  also  any  equipment  errors  in  design  and 
construction,  including  the  effect  of  quantization  of  the 
signal  input.  It  is  not  possible  to  separate  human  errors 
and  equipment  errors,  but  examination  of  the  records 
indicates  that  human  errors  are  almost  certainly  the  larger 


CORRECTIONS  TO  BE  MADE: 
PERISCOPE  ZERO.  ADD +11'. 


-.5  0 +.5 

BEARING  CALIBRATION  <”) 


+ 1.5 


Figure  18.  Distribution  of  errors  in  figure  17.  Correc- 
tions were  made  for  parallax  and  velocity  of  sound. 
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The  inherent  bearing  precision  of  DIMUS  is  believed  to  be 
greater  than  the  results  of  this  test  would  indicate.  However, 
a test  procedure  eliminating  or  greatly  reducing  human 
judgment  would  be  required  to  demonstrate  this. 

The  delay  line  tests  which  were  made  at  the  same  time  as 
the  DIMUS  tests  gave  practically  the  same  results,  namely 
the  same  zero  point  correction  and  standard  deviation. 

Reliability  of  this  DIMUS  system  has  hinged  largely  on  the 
reliability  of  the  shift- register  boards  in  which,  up  to  the 
present,  a fair  number  of  failures  have  occurred,  almost 
solely  due  to  poor  solder  connections.  A preventive  main- 
tenance program  is  planned. 


CONCLUSIONS 


Tests  indicate  that  the  DIMUS  system,  using  five  samples 
per  cycle,  has  satisfactory  beam- forming  capability  and 
its  bearing  precision  is  equal  to  that  obtainable  with  other 
methods,  such  as  delay  lines. 

The  system  may  be  used  in  place  of  delay  lines  in  future 
Lorad  designs,  provided  that  the  simpler,  passive  phase 
compensator  cannot  be  employed. 
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